Navigation between quantum states by quantum mirrors 
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We introduce a technique which allows one to connect any two arbitrary (pure or mixed) su- 
perposition states of an iV-state quantum system. The proposed solution to this inverse quantum 
mechanical problem is analytical, exact and very compact. The technique uses standard and gen- 
eralized quantum Householder reflections (QHR) [Ivanov et al, Phys. Rev. A 74, 022323 (2006)], 
which require external pulses of precise areas and frequencies. We show that any two pure states 
can be linked by two standard QHRs, or by only one generalized QHR. The transfer between any 
two mixed states with the same dynamic invariants (e.g., the same eigenvalues of the density matrix 
p) requires in general N QHRs. Moreover, we propose recipes for synthesis of arbitrary preselected 
mixed states, starting in a single basis state and using a combination of QHRs and incoherent 
processes (pure dephasing or spontaneous emission). 

PACS numbers: 32.80.Bx; 33.80.Be; 03.67.Lx; 03.67.Mn 



I. INTRODUCTION 

Quantum state engineering in atoms and molecules 
traditionally uses three basic techniques for transfer of 
population, complete or partial, from one bound energy 
state to another, single or superposition state: resonant 
pulses of precise areas (e.g. it pulses in a two-state sys- 
tem or generalized 7r pulses for multiple states) adi- 
abatic passage using one or more level crossings 0, or 
stimulated Raman adiabatic passage (STIRAP) and its 
extensions Q . All these techniques require the system to 
be initially in a single energy state; such a state can be 
easily prepared experimentally, e.g. by optical pumping. 
Some of these techniques are "tuned" to a specific initial 
condition: for example, STIRAP requires a counterintu- 
itive pulse sequence to transfer population from state 1 
to 3 in a 1-2-3 linkage, but it is largely irrelevant if the 
system starts in states 2 or 3 (with some exceptions for 
state 3) Q. In other words, STIRAP is (very) useful 
in producing only one column (the first) of the unitary 
propagator. Similar conclusions apply, to a large extent, 
also to the other two techniques using pulse areas and 
level crossings. 

These traditional techniques resolve only a small (al- 
though important) part of the general problem of quan- 
tum state engineering: given the initial and final states of 
an A-state system, find a physical set of operations that 
connect them. This problem requires the construction of 
the entire propagator, not just a single column or row. 

In this paper we introduce a technique for full quantum 
state engineering, which produces in a systematic manner 
a propagator that can connect any two preselected super- 
position states of an A-state quantum system, represent- 
ing a qunit in quantum information [J|. The two states 
can be pure as well as mixed, and the latter may have the 
same or different sets of dynamic invariants (constants of 
motion). The solution consists of two steps: first, find 
a propagator that connects the two states, and second, 
find a physical realization of this propagator. 



The first part is the mathematical solution of this in- 
verse problem in quantum mechanics, and the solution 
is different for three types of problems: (i) pure-to-pure 
states; (ii) mixed-to-mixed states with the same invari- 
ants; (iii) mixed-to-mixed states with different invari- 
ants. The case (iii), for instance, contains the impor- 
tant problem of engineering an arbitrary presected mixed 
state and we pay special attention to it. In this lat- 
ter respect our exact analytic results are alternative to 
the (approximate) numeric optimization procedure pro- 
posed by Karpati et al [5|; moreover, our approach allows 
one to engineer any preselected mixed state, whereas the 
method of Karpati et al [e| can only produce a class of 
mixed states. 

The second part of the solution is the physical real- 
ization of the respective propagator. For this we use the 
recently introduced physical implementation of the quan- 
tum Householder reflection (QHR) || 0] and we show 
that QHR is a very powerful tool for quantum state en- 
gineering. Remarkably, in case (i) only a single QHR is 
needed to connect two pure states. In case (ii), a general 
U(A) propagator is necessary in the general case, which 
requires A QHRs. In case (iii), some sort of incoherent 
process is required in order to equalize the different dy- 
namic invariants of the initial and final mixed states, and 
the remaining coherent U(A) part is realized by QHRs. 
We describe the use of two such incoherent processes: 
pure dephasing and spontaneous emission. 

The Householder reflection Q is a powerful and nu- 
merically very robust unitary transformation, which has 
many applications in classical data analysis, e.g., in solv- 
ing systems of linear algebraic equations, finding eigen- 
values of high-dimensional matrices, least-square opti- 
mization, QR decomposition, etc. [9(. In its quantum 
mechanical implementation p, @] it consists of a single 
interaction step involving N simultaneous pulsed fields 
of precise areas and detunings in an A-pod linkage pat- 
tern, wherein the A states of our system are coupled to 
each other via an ancillary excited state, as displayed 
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the normalization condition and the unimportant global 
phase accounted for). 

The generalized QHR is defined as 



M(v;cp) = I + 



1 



(2) 



where y> is an arbitrary phase. The standard QHR flU) 
is a special case of the generalized QHR ([2) for = 
7r: M(u;7r) = M(w). The generalized QHR is unitary, 
M.(v; (p)^ 1 = M(u; (/j)' = M(u; — </?), and its determinant 
is detM = e iv . 



FIG. 1: Physical realization of the quantum Householder re- 
flection: N degenerate (in RWA sense) ground states, forming 
the qunit, coherently coupled via a common excited state by 
pulsed external fields of the same time dependence and the 
same detuning, but possibly different amplitudes and phases. 



in Fig. [TJ We use two types of QHRs: standard and 
generalized; the latter involves an additional phase fac- 
tor. The standard QHR can operate on or off resonance, 
whereas the generalized QHR requires specific detunings. 
Any unitary matrix can be decomposed into (and there- 
fore, synthesized by) N — I standard QHRs and a phase 
gate, or into N generalized QHRs, without a phase gate; 
hence only N physical operations are needed, which al- 
lows one to greatly reduce the number of physical steps, 
from 0(N 2 ) in existing U(2) realizations [IfJ| to only 
O(N) with QHRs. 

This paper is organized as follows. In Sec. |TT]we re- 
view the standard and generalized QHR gates and their 
physical implementations. In Sec. IIIII we show how two 
pure states can be connected by means of standard and 
generalized QHRs. In Sec. |lV]we construct the propaga- 
tor connecting two arbitrary mixed states with the same 
dynamic invariants. Engineering of an arbitrary prese- 
lected mixed qunit state is presented in Sec. El The 
conclusions are summarized in Sec. IVII 



II. THE TOOL: QUANTUM HOUSEHOLDER 
REFLECTION 



Definition 



The standard QHR is defined as 



M(v) =1-2 



v {v 



(1) 



where I is the identity operator and \v) is an N- 
dimensional normalized complex column- vector. The 
QHR ([TJ) is both hermitian and unitary, M = M f = 
M _1 , which means that M is involutary, M 2 = I. In ad- 
dition, detM = — 1. For real \v) the Householder trans- 
formation |T]) has a simple geometric interpretation: re- 
flection with respect to an (N— l)-dimensional plane with 
a normal vector \v). In general, the vector \v) is complex 
and it is characterized by 2N — 2 real parameters (with 



B. Physical implementation 

We have shown recently @, 0] that the standard and 
generalized QHR operators can be realized physically in 
a coherently coupled iV-pod system shown in Fig. [TJ 
The N degenerate [in the rotating-wave approximation 
(RWA) sense fl] ground states \n) (n = 1,2,..., N), 
which represent the qunit, are coupled coherently by N 
external fields to an ancillary, excited state \e) = \N + 1) 
The excited state \e) can generally be off resonance 
by a detuning A (t) [6| , which must be the same for all 
fields. The Rabi frequencies fii(t), . . . , f2jv(i) of the cou- 
plings between the ground states and the excited state 
have the same pulse-shaped time dependence / (t) , but 
possibly different phases (3 n and amplitudes \ n , 



n n {t)= X nf(t)e^ (n = 1,2,..., AT). 
The qunit +ancilla RWA Hamiltonian reads 



(3) 



HM-! 







n x (t) 
n 2 (t) 



Q N (t) 

n* N (t) 2A(t) 



(4) 



The exact solution to the Schrodinger equation for the 
propagator XJ(t), 



ih—XJ(t)=H(t)V(t). 

at 



(5) 



can be found in [fj . 

The standard QHR M(i>) is realized on exact resonance 
(A = 0), for any pulse shape f(t), and for root-mean- 
square (rms) pulse area 



yl = 2 (2fc + 1) 7T (fc = 0,l,2, 



where 



A = 



ti(t)dt, 



(6) 



(7) 



with n(t) 



1/2 



Then the transition 



probabilities to the ancilla state vanish and the prop- 
agator within the qunit space is given exactly by the 
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standard QHR M(v) (PTJ) . The components of the N- 
dimensional normalized complex vector \v) are the Rabi 
frequencies, with the accompanying phases Q, 



with the desired areas. Of course, the use of femtosecond 
pulses offers another advantage: decoherence is irrelevant 
on such time scales. 



X 



1/2 



(8) 



where X — (j2n=iXnJ ■ Hence the qunit propagator 
represents a physical realization of the standard QHR in 
a single interaction step. Any QHR vector (|8|) can be 
produced on demand by appropriately selecting the peak 
couplings \n and the phases (3 n of the external fields. 

The generalized QHR M(u; (p) can be realized in the 
same TV-pod system, but for specific detunings off reso- 
nance. Again the transition probabilities to the ancilla 
state must vanish; the corresponding rms pulse areas ([7]) 
in general depend on the pulse shape and differ from the 
resonance values ©. The propagator within the qunit 
space is the generalized QHR ©, wherein the phase (p 
depends on the interaction parameters. Although the 
parameters (i.e. the rms area and the detuning) of any 
needed generalized QHR can be found numerically for 
essentially any pulse shape, it is very convenient to use a 
hyperbolic-secant pulse shape, for which there is a simple 
exact analytic solution: the Rosen-Zener model [HI], 



/(*) 
A(t) 



sech (t/T) . 
Ac 



(9a) 
(9b) 



For this pulse shape, the rms area Q is A = nxT. A 
generalized QHR transformation M(v; ip) is realized 
when the interactions satisfy again Eq. ((8]), and the pulse 
area and the detuning obey [fl 0] 



A = 2nl (1 = 1,2,...). 



= 2argJJ[A T + i(2fc+l)]. 



(10a) 
(10b) 



fe=0 



For any given if, there are / values of A , which satisfy 
Eq. (|10b[) 6]. This is also the case for ip — tt, i.e. for the 
standard QHR, for which one of the solutions is Ao = 0. 
Hence the standard QHR M(w) can be realized both on 
and off resonance, whereas the generalized QHR M(w; ip) 
can only be realized for nonzero Ao. The advantage of 
tuning off resonance is the lower transient population in 
the ancilla excited state, which would reduce the popula- 
tion losses if the lifetime of this state is short compared 
to the interaction duration. 

This implementation is particularly suited for a qutrit 
(N = 3) formed of the magnetic sublevels of an atomic 
level with angular momentum J = 1; then the ancilla 
state should be a J = state. The three pulsed fields 
can be delivered from the same laser by using beam split- 
ters and polarizers, which would authomatically ensure 
that all of them have the same detuning and pulse shape. 
Moreover, with femtosecond pulses it would be possible 
to use pulse shapers [12[ , which can easily deliver pulses 



C. Householder decomposition of a U(iV) 
propagator 

The standard QHR M(w) and the generalized QHR 
M.(v;ip) can be used for U(N) decomposition Q- Any 
iV-dimensional unitary matrix U (U _1 = TJt) can be 
expressed as a product of N — 1 standard QHRs M (v n ) 
(n = 1, 2, . . . N — 1) and a phase gate, 

N 

$>(ct> 1 ,...,cj> N ) = J2 ^" l») H = diag {e* , . . . , } , 

n=l 

(11) 

as 

U = M {vi) M (v 2 ) ...M (vn-i) $ fa, ■ ■ ■ , <Pn) , 

(12) 

or as a product of N generalized QHRs, 

U = M (vi;tpi) M (v 2 ; <p 2 ) . . .M (v N ; <p N ) . (13) 



III. TRANSITION BETWEEN PURE STATES 

The designed recipe for constructing a general U(7V) 
transformation makes it possible to solve the important 
quantum mechanical problem of transfering an iV-state 
quantum system from one arbitrary preselected initial 
superposition state to another such state, i.e. the inverse 
problem of quantum-state engineering. The cases of pure 
and mixed states require separate analyses. 



A. Transition by standard QHR 

1. General case 

A pure qunit state is described by a state vector = 
Yln=l c n\ n )i where the vectors \n) represent the qunit 
basis states, and c„ is the complex-valued probability 
amplitude of state \n). Given the preselected initial state 
and the final state Y$> f) of the qunit, we wish to find 
a propagator U, such that 



|*/>=U|*< 



(14) 



We shall show that one of the possible solutions of Eq. 
I4l) reads 



U = Mfu/jDMf!), 



(15) 



where M(wi) and M(w/) are standard QHRs. Here D 
is an A-dimensional unitary matrix, which, when acting 
upon a single qunit basis state only shifts its phase, 



Din) 



e**» In) 



(16) 
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For example, D can be an arbitrary TV-dimensional phase 
gate (fTTj) . Alternatively, D can be an arbitrary TV- 
dimensional generalized QHR M (v; ip) with vector \v) 
orthogonal to the qunit state |n), (v \n) = 0. Finally D 
can be the identity D = I. 

In order to prove Eq. (fT5]) we first define the vector 



(17) 



P 

? 3 



V2[l-Re((* Q |n) e^-»)]' 



where \n) is an arbitrarily chosen basis qunit state, (p an = 
arg[* a ] n and a = i,f. The QHR M.(v in ) acting upon 
|*j) reflects it onto the single qunit state \n), 



M(v m ) = e* 



(18) 



The action of D upon \n) only shifts its phase, see Eq. 
(fTB"]) . The action of M(u/ n ) upon \n) reflects this vector 
onto the final state, 

M(v fn ) \n) = e"**'- |#,) . (19) 

Equations (JTSJl, JTBJ) and (TTH) imply that 

M(v fn )T>M(v in ) |*0 = e <(*'«--*'/-+*») 1^) , (20) 

which, up to an unimportant phase, proves Eq. (|15p . 

The arbitrariness in the choice of the unitary matrix 
D, and the intermediate basis state \n), means that the 
solution (|15p for U is not unique. However, what is im- 
portant is that it always exists. In fact the availability 
of multiple solutions offers some flexibility for a physi- 
cal implementation. In particular we can always choose 
D = I; then the physical realization of the propagator U 
requires only two standard QHRs. 

2. Special cases 

In several important special cases only a single stan- 
dard QHR is needed for pure-to-pure transition. 

1. If the qunit is in a single initial basis state 1*0 = 
|n) then, as follows from Eq. (fT9|) . only one standard 
QHR M(t>/ n ) is sufficient to transfer it into an arbitrary 
superposition state I*/), with |u/ n ) given by Eq. ((T7|) . 

2. Likewise, an arbitrary inital superposition state |*0 
can be linked to any single final state |*/) = \n) by only 
one standard QHR M.(v in ), with |u in ) given by Eq. ([17]) . 

3. If |*0 and |*/) are orthogonal ((*/|*0 = 0), 
then again only a single standard QHR M(u), with \v) — 

(I*/) — 1^))' i s sufficient to connect them. 

4. If |*) and |*^) are superpositional states with 
real coefficients, then again a single standard QHR M(w), 

with \v) = (|*/) - I*,)) /V 2 (!- W7Wdj, is sufficient 
to link them. 











JJ a> 






Time (units of T) 

FIG. 2: (color online) Time evolution of the pulsed fields 
(top) and the populations (bottom) for the transition (I23[) in 
a qutrit. We have assumed sech pulses and rms pulse area 
A — 4% (xT = 4). The individual couplings Xn (n = 1, 2, 3) 
are given by the components of the QHR vector (|24|) . each 
multiplied by \. The detuning is AoT = 1.732 (which gives 
ip — 7r). The thick curve is the state mismatch (|25[) . 



B. Transition by generalized QHR 

Generalized QHR is ideally suited for a pure-to-pure 
transition because, as it is easily seen, only one general- 
ized QHR is sufficient to reflect state |*0 onto |*/), 



where 



\v) 



U = M(y;<p), 



\*f)-m 

2arg(l- <* / |* i )) + 7 r. 



(21) 

(22a) 
(22b) 



In comparison with (|15jl the solution pip is unique; there 
is no arbitrariness in the choice of the QHR vector |u) (up 
to an unimportant global phase) and the phase ip. 



C. Examples 

We consider a qutrit (N = 3), for which the QHR im- 
plementation is particularly suitable. As a first example, 
the transition from a single state to a superposition state, 
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|1) + |2) + |3) 



I** 



(23) 
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FIG. 3: (color online) Time evolution of the pulsed fields 
(top) and the populations (bottom) for the transition (|26[) 
in a qutrit. We have assumed sech pulse shapes and rms 
pulse area A = 2n (\T — 2). The individual couplings Xn 
(n = 1,2,3) are given by the components of the standard 
QHRs (J27J, each multiplied by \- The detuning is AT = 0. 
The thick curve shows the state mismatch (12511, 



FIG. 4: (color online) Time evolution of the pulsed fields 
(top) and the populations (bottom) for the transition (|26[) 
in a qutrit. We have assumed sech pulse shapes and rms 
pulse area A = 2tt (\T — 2). The individual couplings x™ 
(n = 1, 2, 3) are given by the components of the generalized 
QHR HD, each multiplied by x- The detuning is AT = 0.791, 
which produces the desired phase ip = 0.5747T. The thick 
curve shows the state mismatch (I25|l . 



is performed by a single QHR M(t>), with 



V3- 1,-1,-1 




(24) 



Figure [2] shows the corresponding time evolutions of the 
populations and the state mismatch D. The latter is 
defined as the distance between the qutrit state vector 
!*(£)) and the target state I*/), 



D(t) 



^2mn \Pmn(t) 



P f I 



E, 



(25) 



where p m n are the elements of the qutrit density matrix 
p. This definition of D applies to pure and mixed states 
as well. The behavior of D allows us to verify that not 
only the populations but also the phases of the probabil- 
ity amplitudes of the target state \^ f) are produced by 
the QHR. Indeed, as time progresses, D approaches zero, 
which implies that | *&(£)) aligns with l^/). 

In another example, we transfer a two-state superpo- 
sition to a three-state superposition, 



|1) + |3) 
v/2 



ID 



„»/3 



|2) 



Jir/7 



|3) 



V3 



(26) 



by two standard QHRs, V = M(v f )M(vi), with 



= [-0.383, 0,0.924f 



-0.460, 0.628e i7r/3 , 0.628e l7r/7 



(27a) 
(27b) 



or by one generalized QHR, U = M(t>; cp), with 



|u) = [0.194e a2137ri ,0.863e" 



-0.4547ri 



0.467e 



-0.0837ril 



(28) 

and ip — 0.5747T. Figure [3] shows the time evolution 
of the populations and the state mismatch ([25]) for a 
standard-QHR implementation, and Fig. 3] for gener- 
alized QHR. In both cases, the mismatch D vaniches, 
indicating the creation of the desired superposition (|26[) . 
The generalizcd-QHR implementation is clearly superior 
because it creates the target state in a single step. 

In conclusion of this section, we have demonstrated 
that any two pure superposition qunit states can be con- 
nected by just a single generalized QHR, or by two stan- 
dard QHRs. This suggests that QHR, and particularly 
the generalized version, is the most convenient and ef- 
ficient tool for pure-to-pure state navigation in Hilbert 
space. We now turn our attention to the problem of con- 
necting two arbitrary mixed states. 
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IV. COHERENT NAVIGATION BETWEEN 
MIXED STATES 

A mixed qunit state can be described by its density 
matrix p, whose spectral decomposition reads 



N 



(29) 



The eigenvalues r n of p satisfy Y^,n=i r ™ = -*-> ana - 
are the orthonormalized ((^ |^ n ) = Skn) complex eigen- 
vectors of p. The density matrix is hermitian; hence it 
can be parameterized by N 2 — 1 real parameters. 

A hermitian Hamiltonian induces unitary evolution be- 
tween an initial mixed state pi and a final state pf, 



Pf = U Pl Ut. 



(30) 



A unitary evolution does not change the eigenvalues 
{ r «}n=i' wn i cn are therefore dynamic invariants, which is 
easily seen from Eq. (|3"Tj)) (as an equivalent set of dynamic 
invariants one can take the set {Trp™}^ =1 ). Therefore, a 
unitary propagator U can only connect mixed states with 
the same set of invariants {r n } v In order to connect 
mixed states with different invariants we need an incoher- 
ent process; we shall return to this problem in the next 
section. Here we shall find the solution to the problem 
of linking two mixed states with the same invariants. 

Because the eigenvalues {r n } , of pi and pj are the 
same, we should have 



R-PjRj = Rtp/R 



fPf^f - Po, 



(31) 



where the unitary matrices R; and Ry diagonalize re- 
spectively pi and pf, and po =diag{ri, r^, ■ ■ ■ , r^}- By 
replacing Eq. dSTJ) into Eq. $5$ we find 



po = T>p T)\ 

r}ur, 



D 



(32a) 
(32b) 



Because D is a unitary matrix we find poT> = Dpo. Since 
po is diagonal, D must be diagonal too. There are no 
other restrictions on D; hence D is an arbitrary diagonal 
matrix. It follows from Eq. (|32b[) that the solution for 
the unitary propagator in Eq. (f3"0"]) is given by 



U = R/DR 



(33) 



The propagator (|33[) is not unique; it depends on the 
choice of the diagonal matrix D. In particular, we can 
always choose D = I. Hence, the transfer between two 
mixed states requires a general U(iV) propagator. The 
latter can be expressed as a product of N — 1 standard 
QHRs M(u„) (n = 1, 2, . . . , N — 1) and a phase gate 
<fr ((j>i, cf>2, ■•■ , 4>n), Eq. (fT2"j). or by N generalized QHRs 
M{v n ;v n ) {n=l,2,...,N), Eq. ffi 



We take as an example a qutrit, with the arbitrarily 
chosen initial and final density matrices 



Pi = 



Pf = 



0.490 
O.H5e°- 789 « 
0.158e-°- 107 ™ 



0.115e~ - 789 ™ 

0.336 
O.OlSe - 675 ™ 



0.158e 1077 ™ 
0.018e-°- 675,H 
0.175 



(34a) 



0.298 0.022e - 6897 " 0.033e - 319 ™ 
0.022e-°- 689 " 0.180 0.177e 09097 " 
0.033e-°- 319 ™ 0.177e-°- 909 ™ 0.523 



(. 34b ) 

These density matrices can be connected by the unitary 
propagator (f3"3"|) with D = I: U = R^Rj. The latter can 
be expressed as a product of two standard QHRs and one 
phase gate U = M(i>i)M(u 2 )*, with 



K) = [0.612e a532 ™ 0.091e a211 ™ 0.785e a690 ™ 



\»2) 



= [0,0.533 
diag{e 



e -o.i8i7ri 0>846e 



-0.4687ri „0.8197ri 



0.8597ril 



-0.350™ 



}• 



|(§5a) 
(35b) 
(35c) 



or by three generalized QHRs, 
M.{v 1 ;<p 1 )M(v 2 ;>p2)M-(v 3 ;ip 3 ), with 



U 



[0.721e° 



\v 2 ) = [0,0.813e° 



4697ri 



0.080e- - 2097rl ,0.689e°- 270 ^3 T 6a) 

0.582e' a26l7rl ] T , (36b) 

v 3 ) = [0,0, 1] T , (36c) 
Vl = -0.8417T, ip 2 = 0.969tt, ip 3 = -0.128tt. (36d) 

Figure [5] shows the respective time evolution of 
the populations and the state mismatch (|2"5)) for the 
generalized-QHR realization (|36|) . The first QHR 
M(t>3;iy93) does not cause population changes because it 
is in fact a phase gate. As time progresses, the mis- 
match decreases and the target density matrix (134b)) is 
approached. 



V. SYNTHESIS OF ARBITRARY 
PRESELECTED MIXED STATES 

As it was shown in the previous sections, by applying 
one or more QHRs one can connect any two arbitrary 
pure states, or two arbitrary mixed states with the same 
dynamic invariants {f* n } n=1 . Mixed states with different 
invariants cannot be connected by coherent hermitian 
evolution because these invarants are constants of mo- 
tion. Hence in order to connect mixed states with differ- 
ent invariants we need a mechanism with non-hermitian 
dynamics, which can alter the dynamic invariants. 

In this section we shall describe two techniques for en- 
gineering an arbitrary mixed state, starting from a single 
pure state. This is the most interesting special case of 
the general problem of connecting two arbitrary mixed 
states, because the initial state can be prepared rou- 
tinely by optical pumping. Moreover, the general mixed- 
to-mixed problem can be reduced to the single-to-mixed 
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FIG. 5: (color online) Time evolution of the pulsed fields 
(top), the populations and the state mismatch (bottom) for 
the transition between states (|34ap and (|34b|) in a qutrit. We 
have assumed sech pulse shapes and rms pulse area A — 2n 
(\T — 2). The individual couplings \n (n = 1,2, 3) are given 
by the components of the generalized QHR (|36|l . each multi- 
plied by x- The detunings are AiT = -0.255, A 2 T = 0.0 49, 
and A 3 T = -4.918, which produce the QHR phases (|36d[l . 



problem by optically pumping the initial mixed state into 
a single state. 

The two techniques use a combination of coherent 
and incoherent evolutions. The coherent evolution uses 
QHRs, whereas the incoherent non-hermitian evolution is 
induced either by pure dephasing or spontaneous emis- 
sion. We shall consider the two techniques separately. 



the density matrix will be diagonal, with the eigen- 
values r n of pf on the diagonal, which implies that 
it will have the same dynamic invariants as, pf. 

• The third step is to connect this intermediate state 
to the desired state pf by a sequence of QHRs, as 
explained in the previous Sec. IIVI 

In summary, we need three steps: a single QHR, a 
dephasing process, and a sequence of N QHRs. Figure 
[5] shows the evolution of the populations and the state 
mismatch (I25[) during the engineering of the mixed state 
(|34b[) by the dephasing technique. The first step is the 
single QHR M(v), with QHR vector 



|v) = [-0.336, 0.816, 0.471] , 

which transfers the single initial state |1) 
superposition state 



Pi 



0.6 V0T8 VO06 
V0T8 0.3 

Vom Vom o.i 



(37) 

to the pure 
(38) 



The second step is the pure dephasing process, which 
nullifies all coherences and leaves the density matrix in a 
diagonal form, 



p 2 = diag{0.6,0.3,0.1}. 



(39) 



The third step is a sequence of two generalized QHRs, 
which transfer p2 into the desired final density matrix 
Pf, Eq. (|34b|) . The QHR components read 



K) = [0.689e a4547rl ,0.280e - 436 ",0.668e-°- 477;:; 



\v 2 ) = [0,0.793e°- 7407rl ,0.609e l 
if i = 0.950tt, ip 2 = -0.7607T 



0.0257rq 



(40b) 
(40c) 



A. Using dephasing 

We assume that the qunit is initially in the single qunit 
state pi = \i) and we wish to transform the system to 
an arbitrary mixed state pf. Let us denote the eigenval- 
ues of pf by r n (n = 1, 2, . . . , N). We proceed as follows. 

• First, using the prescription from Sec. IIII1 we apply 
a single QHR to transfer state \i) to a pure superpo- 
sition state, in which the populations are equal to 
the eigenvalues of pf\ p nn = r n (n = 1, 2, . . . , N). 
The phases of this superposition are irrelevant. 

• In the second step we switch the dephasing on 
and let all coherences decay to zero. This can be 
done, for example, by using phase-fluctuating far- 
off-resonance laser fields. In the end of this process, 



B. Using spontaneous emission 

In the method, which uses spontaneous emission, we 
start again in a single qunit state pi = \i) {i\, and the 
target is the arbitrary mixed state pf. The procedure 
now consists of only two steps: incoherent and coherent. 
It is particularly well suited for a qutrit, which we shall 
describe, although it is readily extended to more states. 
This method requires a closed qunit-ancilla transition; if 
the ancilla state can decay to other levels then the fidelity 
will be reduced accordingly. 

It is possible here to apply directly the incoherent step, 
which produces a density matrix with the desired final 
dynamic invariants, without the need to prepare first a 
coherent qunit superposition, as in the dephasing method 
above. The idea is to use laser-induced spontaneous emis- 
sion from the ancilla excited state to prepare a completely 



8 



rt 2 



.2 1 



M(v;0 deph. M(v 2 ;«,) M(v ;p) 





















i 


ill 












-10 10 20 
Time (units of 7) 



40 



FIG. 6: (color online) Time evolution of the pulsed fields 
(top), the populations and the state mismatch (I25|) (bottom) 
for mixed state engineering in a qutrit. The qutrit starts in 
state |1) and the target final state is given by Eq. ()34b[l . 
We have assumed sech pulse shapes and rms pulse area A = 
2n {\T — 2). The individual couplings Xn (n = 1,2,3) are 
given by the components of the generalized QHR (|40p . each 
multiplied by \- Tne detunings are AiT = 0.072 and A 2 T = 
—0.396, which produce the desired QHR phases (I40c|l . The 
dephasing rate is T = 2/T. 

incoherent superposition of the qunit states with popu- 
lations p nn equal to the eigenvalues r n of pj, 



(41) 



For this we apply a sequence of appropriately chosen laser 
pulses from the qunit states to the excited state, which 
decays back to the qunit states and redistridutes the pop- 
ulation among them. 

There are various scenarios possible, which can pro- 
duce the desired incoherent qunit superposition. Here we 
describe a scenario which looks particularly simple and 
easy to implement for the qutrit formed of the magnetic 
sublevels M = —1,0, 1 of a J = 1 level and an ancilla 
excited level with J = (this implies also equal sponta- 
neous decay branch ratios from the J = level to the 
M sublevels of the qutrit). For definiteness, and with- 
out loss of generality, we assume that the eigenvalues of 
Pf are ordered as t\ ^ r 2 ^ r%. We need three pulses: a 
short pulse from state |1), a long pulse from state |3) and 
again a short pulse from state |1) (here short and long 



are related to the lifetime of the excited state). 

The short pulse from the initially populated state |1), 
with excitation probability p\ , transfers population p\ to 
the excited state, 1 /3 of which decays back to each of the 
qutrit states. The ensuing density matrix reads 



j | 2 1 1 

Pi = diag <J 1 - -pi, - Pl , -pi 



(42) 



We then apply a sufficiently long pulse from state |3), so 
that its population is completely depleted and distributed 
among states |1) and \2). The resulting density matrix is 



I':: = dm-- { 1 - ipi, ipi,0 



(43) 



We now apply again a short pulse from state |1), with 
a different probability p 2 , and then wait for spontaneous 
emission from the excited state. The result is 



P3 



diag 
1 



1 



1 



1 



7P2 1 



3 

1 

' 3 



P2 



Mi 



.(44) 



It is easy to show that in order to create the mixed state 
(|4l"j) we should have the probabilities 



Pl = 2 (r 2 - r 3 ) , 
3r 3 

P2 = 



(45a) 
(45b) 



Because we assumed that n ^ r 2 ^ r 3 the probabili- 
ties pi and p2 belong to the interval [0, 1] and are there- 
fore well defined. Such probabilities can be produced by 
resonant pulses with appropriate pulse areas A n . These 
pulses should be short compared to the lifetime of the ex- 
cited state in order to avoid spontaneous emission during 
their action. 

Once we have prepared the mixed qutrit state (|41[) . 
which has the same invariants as p/, we can apply QHRs 
to transfer this state into the desired final state pf, as 
described in Sec. HVl 



VI. CONCLUSIONS 

In this paper we have proposed a technique, which al- 
lows to connect any two quantum superposition states, 
pure or mixed, of an N-state atom. This solution of 
the inverse problem in quantum mechanics contains two 
stages: (i) mathematical derivation of the propagator 
that links the desired initial and final density matrices, 
and (ii) physical realization of this propagator. In the 
most general case of arbitrary mixed states, the imple- 
mentations combine coherent hermitian and incoherent 
non-hermitian interactions induced by pulsed laser fields. 
In general, the propagator is not unique, which reflects 
the multitude of paths between two qunit states; this also 
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allows for some flexibility in the choice of most convenient 
path. 

The physical realization uses an iV-pod configuration 
of N lower states, forming the qunit, and an ancillary 
upper state. It is particularly convenient for a qutrit, 
where the N = 3 states arc the magnetic sublevels of 
a J = 1 level and the ancilla state is a J = level. 
Then only a single tunable laser is needed to provide the 
necessary polarized laser pulses. 

The hermitian part uses a sequence of sets of short 
coherent laser pulses with appropriate pulse areas and 
detunings. For each set, the propagator of the iV-pod 
represents a quantum Householder reflection (QHR). A 
sequence of at most N suitably chosen QHRs can syn- 
thesize any desired unitary propagator. 

We have shown that two arbitrary preselected pure su- 
perposition states can be connected by a single QHR 
only, because the respective propagator has exactly the 
QHR symmetry. Two mixed states, with the same set of 
dynamic invariants, require a general U(iV) transforma- 
tion, which can be realized by at most N QHRs. This is a 
significant improvement over the existing setups involv- 
ing 0(N 2 ) operations, which can be crucial in making 
quantum state engineering and operations with qunits 
experimentally feasible. 

The most general case of two arbitrary mixed states 
with different dynamic invariants requires an incoherent 
step, which equalizes the invariants of the initial den- 



sity matrix to those of the final density matrix. We have 
demonstrated how this can be done by using pure dephas- 
ing or spontaneous decay of the ancillary upper state. 
Once the invariants are equalized, the problem is reduced 
to the one of connecting two mixed states with the same 
invariants, which, as explained above, can be done by at 
most N QHRs. This method has been described for a 
qutrit, but it is easily generalized to an arbitrary qunit. 

The present results can have important applications 
in the storage of quantum information. For example, a 
qubit can encode two continuous parameters: the popu- 
lation ratio of the two qubit states and the relative phase 
of their amplitudes. A qunit in a pure state can encode 
2(N— 1) parameters (N—l populations and N—l relative 
phases), i.e. by using qunits information can be encoded 
in significantly fewer particles than with qubits. More- 
over, a mixed qunit state can encode as many as iV 2 — 1 
real parameters. This may be particularly interesting if 
the number of particles that can be used is restricted, 
e.g., due to decoherence Q. 
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